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We calculate the SIN and SIS tunneling conductances for the spin-fermion model. We argue that 
at strong spin-fermion coupling, relevant to cuprates, both conductances have dip features near the 
threshold frequencies when a tunneling electron begin emitting propagating spin excitation. We 
argue that the resonance spin frequency measured in neutron scattering can be inferred from the 
tunneling data by analyzing the derivatives of SIN and SIS conductances. 
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The electron tunneling experiments are powerful tools 
to study the spectroscopy of superconductors. These ex- 
periments measure the dynamical conductance dl/dV 
through a junction as a function of applied volt- 
age V and temperature ||||. For superconductor- 
insulator-normal metal (SIN) junctions, the measured 
dynamical conductance is proportional to the elec- 
tron density of states (DOS) in a superconductor 
N(u) = -(l/ir)fdk ImG(k,u>) at ui = eV [|. 
For superconductor-insulator-superconductor (SIS) junc- 
tions, the conductance dl/dV oc G(u> — eV), where 
G(u) — L d£lN(uj — fi) 9oiV(0) is proportional to the 
derivative over voltage of the convolution of the two 
DOS g. 

For conventional superconductors, the tunneling ex- 
periments have long been considered as one of the most 
relevant ones for the verification of the phononic mech- 
anism of superconductivity Q. In this communication 
we discuss to which extent the tunneling experiments on 
cuprates may provide the information about the pair- 
ing mechanism in high-T c superconductors. More specif- 
ically, we discuss the implications of the spin-fluctuation 
mechanism of high-temperature superconductivity on the 
forms of SIN and SIS dynamical conductances. 

The spin-fluctuation mechanism implies that the pair- 
ing between electrons is mediated by the exchange of 
their collective spin excitations peaked at or near the an- 
tiferromagnetic momentum Q. This mechanism yields a 
d— wave superconductivity H, and explains |s|j7j a num- 
ber of measured features in superconducting cuprates, in- 
cluding the peak/dip/hump features in the ARPES data 
near (0, it) ||, and the resonance peak below 2A in the 
inelastic neutron scattering data ||. Moreover, in the 
spin-fluctuation scenario, the ARPES and neutron fea- 
tures are related: the peak-dip distance in ARPES equals 
the resonance frequency in the dynamical spin suscepti- 
bility . This relation has been experimentally verified 
in optimally doped and underdoped YBCO and opti- 
mally doped Bi2212 materials Here we argue that 
the resonance spin frequency can also be inferred from 
the tunneling data by analyzing the derivatives of SIN 
and SIS conductances. 

The SIN and SIS tunneling experiments have been 
performed on YBCO and 5z2212 materials @||. At 
low/moderate frequencies, both SIN and SIS conduc- 



tances display a behavior which is generally expected in 
a d— wave superconductor: SIN conductance is linear in 
voltage for small voltages, and has a peak at eV = A 
where A is the maximum value of the d— wave gap [jp, 
while SIS conductance is quadratic in voltage for small 
voltages, and has a near discontinuity at eV — 2 A [0. 
These features have been explained by a weak-coupling 
theory, without specifying the nature of the pairing in- 
teraction pl]| . However, above the peaks, both SIN and 
SIS conductances have extra dip/hump features which 
become visible at around optimal doping, and grow with 
underdoping |l],|2| . We argue that these features are sen- 
sitive to the type of the pairing interaction and can be 
explained in the spin-fluctuation theory. 

As a warm-up for the strong coupling analysis, consider 
first SIN and SIS tunneling in a d— wave superconductor 
in the weak coupling limit. In this limit, the fermionic 
self-energy is neglected, and the superconducting gap 
does not depend on frequency. For simplicity, we con- 
sider a circular Fermi surface for which A^ = A cos 2(f). 

We begin with the SIN tunneling. Integrating 
G(k,u>) = (oj + e k )/(uj 2 - t\ - A|) over e k =v F (k- k F ) 
we obtain 



N{uj) = Re 



2tt Jo ,/uj 2 - A 2 cos 2 (20) 

2 f K(A/u) for uj > A 

{u/A)K{w/A) for uj < A 



(1) 



where K(x) is the elliptic integral. We see that 
N(u>) ~ ui for uj <C A and diverges logarithmically as 
(l/V) ln(8A/|A — uj\) for w« A. At larger frequencies, 
N(u>) gradually decreases to a frequency independent, 
normal state value of the DOS, which we normalized to 
1. The plot of N(uj) is presented in Fig |l|a. 

We now turn to the SIS tunneling. Substituting the 
results for the DOS into G{uj) and integrating over O, we 
obtain the result presented in Fig[j]b. At small uj, G(uj) is 
quadratic in frequency, which is an obvious consequence 
of the fact that the DOS is linear in uj. At uj = 2A, G(u) 
undergoes a finite jump. This discontinuity is related to 
the fact that near 2 A, the integral over the two DOS 
includes the region O ss A where both AT(fi) and N(uj — 
SY) are logarithmically singular, and dnN(Q) diverges as 
1/(0 — A). The singular contribution to G(uj) from this 
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FIG. 1. The behavior of SIN and SIS tunneling conduc- 
tances, dl/dV, in a d— wave BCS superconductor (figures a 
and b, respectively) 



region can be evaluated analytically and yields 



1 



-P 



dx In \x\ 



1 



sign(w - 2A) (2) 



We see that the amount of jump in the SIS conductance 
is a universal number which does not depend on A. 

The results for the SIN and SIS conductances in a 
rf— wave gas agree with earlier studies jy]. In previous 
studies, however, SIS conductance was computed numer- 
ically, and the universality of the amount of the jump at 
2A was not discussed, although it is clearly seen in the 
numerical data. 

We now turn to the main subject of the paper and 
discuss the forms of SIN and SIS conductances for strong 
spin-fcrmion interaction. 

We first show that the features observed in a gas are in 
fact quite general and are present in an arbitrary Fermi 
liquid as long as the impurity scattering is weak. Indeed, 
in an arbitrary d— wave superconductor, 



N(uj) oc Im 



(3) 



where <fi is the angle along the Fermi surface, and F(<j>, ui) 
and E(</>, ui) are the retarded anomalous pairing vertex 
and retarded fermionic self-energy at the Fermi surface 
(the latter includes a bare ui term in the fermionic prop- 
agator). The measured superconducting gap A((j>) is a 
solution of F{<f>, A(0)) = E(0, A(0)). 

In the absence of impurity scattering, imE and ImF in 
a superconductor both vanish at T = up to a frequency 
which for arbitrary strong interaction exceeds A. The 
Kramers-Kronig relation then yields at low frequencies 
ReE((f),u)) oc uj, ReF(cj>,u)) oc (4> - 4> no de) where </w e is 
a position of the node of the d— wave gap. Substituting 
these forms into (|^) and integrating over <fi we obtain 
N(ui) oc ui although the prefactor is different from that 
in a gas. The linear behavior of the DOS in turn gives 
rise to the quadratic behavior of the SIS conductance. 



Similarly, expanding E 2 — F 2 near each of the maxima 
of the gap we obtain E 2 (0,o;) — F 2 (cj),uj) oc (uj — A) + 
B{4> - <?W) 2 , where B > 0. Then 



N(w) oc Re 



Lnin-AI 



(CI — A) 



(4) 



This result implies that the SIN conductance in an ar- 
bitrary Fermi liquid still has a logarithmic singularity at 
eV = A, although its residue depends on the strength 
of the interaction. The logarithmical divergence of the 
DOS causes the discontinuity in the SIS conductance by 
the same reasons as in a Fermi gas. 

In the presence of impurities, the logarithmical singu- 
larity is smeared out, and the DOS acquires a nonzero 
value at zero frequency (at least, in the self-consistent 
T— matrix approximation jl^]). However, for small con- 
centration of impurities, this affects the conductances 
only in narrow frequency regions near singularities while 
away from these regions the behavior is the same as in 
the absence of impurities. 

We now show that a strong spin-fcrmion interaction 
gives rise extra features in the SIS and SIN conductances, 
not present in a gas. The qualitative explanation of 
these features is the following. At strong spin-fermion 
coupling, a d-wave superconductor possesses propagat- 
ing, spin-wave type collective spin excitations near an- 
tiferromagnetic momentum Q and at frequencies below 
2 A. These excitations give rise to a sharp peak in the 
dynamical spin susceptibility at a frequency Cl res < 2 A 
H , and also contribute to the damping of fermions near 
hot spots (points at the Fermi surface separated by Q), 
where the spin-mediated d— wave superconducting gap 
is at maximum. If the voltage for SIN tunneling is such 
that eV = Cl res + A, then an electron which tunnels from 
the normal metal, can emit a spin excitation and fall to 
the bottom of the band (see Fig. ||a) loosing its group 
velocity. This obviously leads to a sharp reduction of the 
current and produce a drop in dl/dV. 

Similar effect holds for SIS tunneling. Here however 
one has to first break an electron pair, which costs the 
energy 2 A. After a pair is broken, one of the electrons 
becomes a quasiparticle in a superconductor and takes an 
energy A, while the other tunnels. If eV — 2A + Cl res , 
the electron which tunnels through a barrier has energy 

~\~ Clyfig , and can emit a spin excitation and fall to the 
bottom of the band. This again produces a sharp drop 
in dl/dV (see Fig. |b). 

In the rest of the paper we consider this effect in more 
detail and make quantitative predictions for the experi- 
ments. Our goal is to compute dl/dV for SIN and SIS 
tunneling for strong spin-fermion interaction. 

The point of departure for our analysis is the set of two 
Eliashberg-type equations for the fermionic self-energy 
E w , and the spin polarization operator Iln- The later is 
related to the dynamical spin susceptibility at the anti- 
ferromagnetic momentum by x l (Qi CI) oc 1 — IIq. The 
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FIG. 2. The schematic diagram for the dip features in 
SIN and SIS tunneling conductances (figures a and b, respec- 
tively). For SIN tunneling, the electron which tunnels from a 
normal metal can emit a propagating spin wave if the voltage 
eV — A + fl re3 where Q r es is the minimum frequency for spin 
excitations. After emitting a spin-wave, the electron falls to 
the bottom of the band which leads to a sharp reduction of 
the current and produces a drop in dl /dV . For SIS tunneling, 
the physics is similar, but one first has to break an electron 
pair, which costs energy 2A. 



same set was used in our earlier analysis of the relation 
between ARPES and neutron data 0. In Matsubara 
frequencies these equations read (S Wm = iE(uJ m )) 



So. — UJ r , 



U °=2 
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This set is a simplification of the full set of Eliashberg 
equations that includes also the equation for the anoma- 
lous vertex F(ui) |fl3|| . As in |Q we assume that near op- 
timal doping, the frequency dependence of F(w) is weak 
at w ~ A relevant to our analysis, and replace F(uj) by a 
frequency independent input parameter F. Other input 
parameters in (|) are the dimensionless coupling constant 
R = g/(vF^ ) and a typical spin fluctuation frequency 
uj s f = (n /4)(vp£ ) /g. They are expressed in terms of 
the effective spin-fermion coupling constant g, the Fermi 
velocity at a hot spot vf, and the magnetic correlation 
length £. By all accounts, at and below optimal doping, 
i? > 1 H], i.e., the system behavior falls into the strong 
coupling regime. 

Strictly speaking, the set (||) is valid near hot spots 
where <f> ~ 4>max- Away from hot spots F(<f) is re- 
duced compared to F. We, however, will demonstrate 
that the new features due to spin-fermion interaction are 
produced solely by fermions from hot regions. 

As in [Q , we consider the solution of (||) for the exper- 
imentally relevant case F 3> Ruj s f when the measured 
superconducting gap A ~ F 2 /(R 2 ui s f) 3> uj s f. In this 
situation, at frequencies ~ A, fermionic excitations in the 



normal state are overdamped due to strong spin-fermion 
interaction. In a superconducting state, the form of the 
spin propagator is modified at low frequencies because of 
the gap opening, and this gives rise to a strong feedback 
from superconductivity on the electron DOS. 

More specifically, we argued in jjj that in a super- 
conductor, Ilo at low frequencies SI <C 2A behaves as 
f2 2 /A, i.e., collective spin excitations are undamped, 
propagating spin waves. This behavior is peculiar to 
a superconductor - in the normal state, the spin ex- 
citations are completely overdamped. The propagat- 
ing excitations give rise to the resonance in x{Qi ^) at 
Sl res ~ (Auj s f) 1/2 < uj sf where ReU(n res ) = 10. This 
resonance accounts for the peak in neutron scattering || . 

The presence of a new magnetic propagating mode 
changes the electronic self-energy for electrons near hot 
spots. In the absence of a propagating mode, an electron 
can decay only if its energy exceeds 3A. Due to reso- 
nance, an electron at a hot spot can emit a spin wave 
already when its energy exceeds A + Sl res . It is essential 
that contrary to a conventional electron-electron scatter- 
ing, this process gives rise to a discontinuity in MS(u) 
at the threshold. Indeed, using the spectral represen- 
tation to transform from Matsubara to real frequencies 
in the first equation in (^|), integrating over momentum 
and neglecting for simplicity unessential q 2 in the spin 
susceptibility, we obtain for u > io t h — A + O res 
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(5) 

We see that ImT,(uj) jumps to a finite value at the thresh- 
old. This discontinuity is peculiar to two dimensions. By 
Kramers-Kronig relation, the discontinuity in ImT, gives 
rise to a logarithmical divergence of i?eS at u> = uj t h- 
This in turn gives rise to a vanishing spectral function 
near hot spots, and accounts for a sharp dip in the 
ARPES data @. 

We now show that the singularity in ReY.(uj) causes 
the singularity in the derivatives over voltages of both 
SIN and SIS conductances d 2 I/dV 2 . Indeed, near a hot 
spot, F((j)) = F(l — X(f) 2 ) where cf) — 4>~4> ma x, and A > 0. 
Then, quite generally, Re~E((f>, ui) oc In \lu — uj t h{(t>)\ where 
uJth{<f) = uj t h + C(f) 2 , and C > 0. Substituting this ex- 
pression into the DOS and differentiating over frequency, 
we obtain after a simple algebra 

F 2 {4>) 



dNjuj) 
duj 



<-d u T,((f),uj)d(f) 

1 &(u>th - u) 
In 3 |w — ujth\ V^th — u 



(7) 



where Q(x) is a step function. We see that dN(uj)/duj has 
a one-sided, square-root singularity at uj = uj t ^ L . Physi- 
cally, this implies that the conductance drops when prop- 
agating electrons start emitting spin excitations. Note 
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FIG. 3. The schematic forms of SIN (a) and SIS (b) tun- 
neling conductances for strong spin-fermion interaction. We 
added small impurity scattering to soften singular features re- 
lated to the sharpness of the Fermi surface (see the text). The 
dip/hump features above the peaks are the strong coupling 
effects not present in a gas. The insets show the derivatives 
of conductances above A for SIN tunneling and 2A for SIS 
tunneling. We argue in the text that these derivatives have 
maxima at voltages eV = f2 = A + f2 res for SIN tunneling 
and eV = Q = 2 A + £l res for SIS tunneling, where n res is the 
resonance spin frequency measured in neutron experiments. 



that the typical 4> which contribute to the singularity are 
small (of order \u t h — UJ \ 1 ^ 2 ), which justifies our assertion 
that the singularity is confined to hot spots. 

The singularity in dN(tu)/dui is likely to give rise to a 
dip in N(lo) at lo > u>±h- The argument here is based 
on the fact that if the angular dependence of L0th(4 > ) 
is weak (i.e., C is small), then T,(io t h) ^S> F(u t h), and 
N(oj t h) reaches its normal state value with infinite neg- 
ative derivative. Obviously then, at u> > u>th, N(u>) goes 
below its value in the normal state and should there- 
fore have a minimum at some lo > ujth- Furthermore, at 
larger frequencies, we solved (JsJ) perturbatively in F(lo) 
and found that N(u>) approaches a normal state value 
from above. This implies that besides a dip, N(lo) should 
also display a hump somewhere above LOth- The behavior 
of the SIN conductance is schematically shown in Fig. |^a. 

Similar results hold for SIS tunneling. The derivative 
of the SIS current, d 2 I/dV 2 ~ dG(uj)/duj, is given by 



8G{uj) 
dui 



d u N(u; - n)d n N(Q,)<m 



(8) 



Evaluating the integral in the same way as for SIN tun- 
neling, we find a square-root singularity at lo = ui^ h = 



d 2 I 
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Q(uth - w) 



dV 2 ./ lo - Q - A In 3 \uj th - fi| ^uo th - w 



(9) 



and lo — f2 
singular. 

Again, it is very plausible that the singularity of the 
derivative causes a dip at a frequency lo > u)^ h , and 
a hump at even larger frequency. We stress, however, 
that at exactly Lo^ h , the SIS conductance has an infinite 
derivative, while the dip occurs at a frequency which is 
somewhat larger than Lo^ h . The behavior of the SIS con- 
ductance is presented in Fig ||. 

Qualitatively, the forms of conductances presented in 
Fig | agree with the SIN and SIS data for YBCO and 
Bi2212 materials |l|,f2||. Moreover, recent SIS tunneling 
data for i?i2212 M indicate that the relative distance 
between the peak and the dip (fl re s/ (2 A) in our theory) 
decreases with underdoping. More data analysis is how- 
ever necessary to quantitatively compare tunneling and 
neutron data. 

To summarize, in this paper we considered the forms 
of SIN and SIS conductances both for noninteracting 
fcrmions, and for fermions which strongly interact with 
their own collective spin degrees of freedom. We ar- 
gue that for strong spin-fermion interaction, the reso- 
nance spin frequency fi res measured in neutron scatter- 
ing can be inferred from the tunneling data by analyzing 
the derivatives of SIN and SIS conductances. We found 
that the derivative of the SIN conductance diverges at 
eV = A + il res while the derivative of the SIS conduc- 
tance diverges at eV = 2 A + £l res , where A is the maxi- 
mum value of the d— wave gap. 
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The singularity comes from the region where Q s» u t h 
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